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The crossover from two- to three-dimensional critical behavior of nearly antiferromagnetic itin- 
erant electrons is studied in a regime where the inter-plane single-particle motion of electrons is 
quantum-mechanically incoherent because of thermal fluctuations. This is a relevant regime for 
very anisotropic materials like the cuprates. The problem is studied within the Two-Particle Self- 
Consistent approach (TPSC), that has been previously shown to give a quantitative description of 
Monte Carlo data for the Hubbard model. It is shown that TPSC belongs to the n — > oo limit of the 
O (n) universality class. However, contrary to the usual approaches, cutoffs appear naturally in the 
microscopic TPSC theory so that parameter-free calculations can be done for Hubbard models with 
arbitrary band structure. A general discussion of universality in the renormalized-classical crossover 
. from d — 2 to d — 3 is also given. 
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PACS numbers: 75.10.Lp, 71.27.+a, 71.10.+X, 74.72.-h 



I. INTRODUCTION 

> 

qq ■ A simple model of itinerant antiferromagnets is provided by electrons on a lattice with short-range repulsion. In 
the low temperature phase, the system is in a Spin Density Wave (SDW) state. In three dimensions, above the 
£N| i transition temperature, the electrons form a so-called nearly antiferromagnetic Fermi liquid. Traditional mean-field 
techniques for studying SDW instabilities of Fermi liquids fail completely in low dimension. In two dimensions for 
example, the Random Phase Approximation (RPA) predicts finite temperature antiferromagnetic transitions while 
this is forbidden by the Mermin- Wagner theorem. Nevertheless, one can study universal critical behavior using various 
forms of renormalization group treatments appropriate either for the strong |l| g| || or the weak-coupling limits [Q || . 
The self-consistent-renormalized approach of Moryia |^| also satisfies the Mermin- Wagner theorem in two dimensions. 
Since cutoff-dependent scales are left undetermined by all these approaches they must be found by other methods. 
*^ . For example, in the strong-coupling limit, the spin-stiffness constant of the non-linear a— model must be determined 
from Monte Carlo simulations. In the weak-coupling case however, Monte Carlo simulations are limited to very small 
O ■ systems, of order 10 x 10 that do not allow one to study much of the critical regime. 

Recently, the Two-Particle Self-Consistent approach was developed to obtain from a microscopic model a quan- 
, titative description of itinerant electrons not only far from phase transitions, but also in the critical regime. It was 
shown j?J that in this approach the Mermin- Wagner theorem is satisfied and that, away from the critical regime, 
the approach gives quantitative agreement with Monte Carlo simulations of the nearest-neighbor and next-nearest 
neighbor Q Hubbard model in two dimensions. Quantitative agreement is also obtained as one enters the narrow 
critical regime accessible in Monte Carlo simulations. The approach is restricted to the one-band Hubbard model 
with on-site interaction, but is valid for arbitrary dispersion relation. The TPSC approach also allows one to study 
the case where the instability of the itinerant electrons is at an incommensurate wave-vector, but in this paper we 
restrict ourselves to the case where the order is at the antiferromagnetic wave vector. The self-consistent-renormalized 
approach of Moryia Q cannot deal with the incommensurate case without a priori information. Even though it has 
the same critical behavior as the TPSC approach it does not allow one to obtain quantitative parameter-free results 
from a microscopic Hamiltonian. 

We first show in full generality that the TPSC approach gives the leading term of the critical behavior in a 1/n 
expansion. In other words, it gives the n — > oo limit of the O (n) model where n = 3 is the physically correct 
(Heisenberg) limit. It will be apparent that there is no arbitrariness in cutoff so that, given a microscopic Hubbard 
model, no parameter is left undetermined. One can go with the same theory from the non-critical to the critical 
regime. 

We then show that the previously studied two-dimensional critical regimes, namely quantum-critical Q and renor- 
malized classical jjj are reproduced here to leading order in 1/n. In the quantum critical regime, one usually dis- 
tinguishes two cases [|J: Model A, where the paramagnetic Fermi surface does not intersect the magnetic Brillouin 
zone, and Model B where it does. This distinction is important in the quantum critical regime because it changes the 
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dynamical critical exponent. In this paper, we also give results on Model C, the case of perfect nesting. In this case, 
the microscopic approach shows that modifications to frequency-independent thermodynamic properties can arise. 
In particular, in the two-dimensional perfect-nesting case the usual exponential dependence of correlation length on 
temperature exp {cst/T) can be modified to be roughly exp {cst/T 3 ) in some temperature region of the renormalized 
classical regime. 

Then we study the renormalized-classical crossover from d — 2 to d = 3 in the highly anisotropic case of weakly 
coupled planes. [^) The general theory of such crossover is given in Appendix D, along with a discussion of universal 
crossover functions. In the main text it is shown that in the highly anisotropic case the crossover can occur in a 
rather unusual regime, namely i|| -C ksT^ <C t± where in (tj_) is the inter (intra) plane hopping integral and T/v is 
the three-dimensional Neel temperature. This regime is unusual because even though one is dealing with an itinerant 
fermion system, the inequality i|| <C kgTpf means that the smallest fermionic Matsubara frequency is larger than the 
dispersion in the parallel direction, making the three-dimensional band structure irrelevant for one-particle properties. 
In the language of Refs. pT| , there is "no coherent band motion" in the parallel direction. Physically, the extent of 
the thermal de Broglie wave packet in the direction perpendicular to the planes is smaller than the distance between 
planes, a situation that does not occur in a usual Fermi liquid since in the isotropic case the inequality ksT <C Ef 
implies that the thermal de Broglie wavelength is much larger than the lattice spacing. Another way of describing 
this i|| <C ksTjy <C t± situation is to say that the itinerant electrons become unstable at the two-particle level while 
their motion in the third direction is still quasi-classical, or quantum incoherent, at the single-particle level because of 
thermal fluctuations. In the more usual situation, coherence at the one-particle level is established before the phase 
transition, namely /csT/v <C in « ii. These two regimes have been extensively discussed in the d = 1 to d = 3 
crossover of Luttinger liquids by Bourbonnais and Caron (!(]] |n[] . 

The above single-particle incoherent regime k^Tj^ <C ij_ is likely to be the relevant one for high-temperature 
superconductors. While the parent insulating compound LaiCuO^ has been extensively studied in the strong coupling 
limit, this type of compound is expected to be in an intermediate-coupling regime. Hence, it is legitimate to approach 
the problem not only from the strong-coupling limit but also from the weak-coupling side, especially with the 
TPSC approach where all cutoffs are determined by the microscopic model. This problem is commented on at the 
end of the paper. More detailed quantitative comparisons with experiment will appear later. 



II. TWO-PARTICLE SELF-CONSISTENT APPROACH 



We start from the Hubbard model, 

H = - ]T tij (ct CT c iCT + c^Cicr) +Uj2 n ^ n il ■ (!) 

<ij>cr i 

In this expression, the operator Ci G destroys an electron of spin a at site i. Its adjoint c\ a creates an electron. The 
symmetric hopping matrix i,j determines the band structure. Double occupation of a site costs an energy U due to 
the screened Coulomb interaction. In the present section, the hopping parameters need not be specified. We work 
in units where fee = 1, and h = 1. As an example that occurs later, the dispersion relation in the d-dimensional 
nearest-neighbor model when the lattice spacing is a is given by 

d 

e k = -2t^2 {cos ha) . (2) 

i=l 

The nearest-neighbor quasi-two dimensional case will be another case of interest later, 

£k = — 2<j_ (cos k x a± + cos k y a±) — 2t\\ cos fc z aii . (3) 



The TPSC approach 0, |13| can be summarized as follows. One approximates spin and charge susceptibilities 
Xspi Xch by RPA-like forms but with two different effective interactions U sp and U c h which are then determined 
self-consistently. Although the susceptibilities have an RPA functional form, the physical properties of the theory 
are very different from RPA because of the self-consistency conditions on U sp and U c h- The necessity to have two 
different effective interactions for spin and for charge is dictated by the Pauli exclusion principle (n^.) = (n a ) which 
implies that both \sp and Xch are related to only one local pair correlation function (n^ni). Indeed, using the 
fluctuation-dissipation theorem in Matsubara formalism we have the exact sum rules, 

(nf) + (nj) + 2(n T nx) - n 2 = ± £ Xch {q) (4) 
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and 



(n*) + (n 2 l )-2(n 1 n L ) = ^}_ / Xsp{q) (5) 

9 

where (3 = 1/T, n = (n-f) + (nj.), g= (q, iq n ) with q the wave vectors of an N site lattice, and with iq n = 2ninT the 
bosonic Matsubara frequencies. The Pauli principle (n 2 a ) = (n a ) applied to the left-hand side of both equations with 
our RPA-like forms for Xsp, Xch on the right-hand side lead to 

n + 2(n T n x ) - n = — ^ : , i „ pr, (6) 

pA ^ 1 + ^UchXoW 

9 

n - 2(n r n l ) = — ^ rr T Fv ( 7 ) 

pN ^ 1- ^U S pXo{q) 

9 

with Xo(<?) the susceptibility for non-interacting electrons. 

If (n-fn^) is known, U sp and C/c^ are determined from the above equations. This key quantity (n^ni) can be obtained 
from Monte Carlo simulations or by other means. However, it may be also be obtained self-consistently by adding 
to the above set of equations the relation 

U sp = g n (0)U ; ^( ) = ^1^!_ ( 8 ) 

Eqs.(0) and (gj) define a set of self-consistent equations for U sp that involve only two-particle quantities. We call 
this ap pro ach Two-Particle Self-Consistent to contrast it with other conserving approximations like Hartree-Fock or 
FLEX p^l that are self-consistent at the one-particle level, but not at the two-particle level. The above procedure Q] 
reproduces both Kanamori-Brucckner screening as well as the effect of Mcrmin- Wagner thermal fluctuations, giving 
a phase transition only at zero-temperature in two dimensions, as discussed in the following section. Quantitative 
agreement with Monte Carlo simulations on the nearest-neighbor |7j and next-nearest-neighbor models || is obtained 
for all fillings and temperatures in the weak to intermediate coupling regime U < St. 

We emphasize that deep in the critical regime, the ansatz Eq.(|^) fails in the sense that gfj.(0) eventually reaches 
zero at T = in the nearest-neighbor Hubbard model at half-filling while there is no reason to believe that this really 
happens. The physically appropriate choice in the renormalized classical regime described below, is to keep the value 
of 3|j.(0) fixed at its crossover-temperature value. In the numerical calculations also described below, we are never 
far enough from Tx to have to worry about this. The value of <?tj,(0) is the one that is determined self-consistently. 



III. CRITICAL BEHAVIOR OF THE TPSC APPROACH IN ARBITRARY DIMENSION 

In this section we discuss the critical behavior of the TPSC approach in arbitrary dimension for hypercubic systems. 
It is convenient to set the lattice spacing to unity. 

As one approaches a phase transition, one enters the renormalized classical regime, [0 where classical thermal 
fluctuations dominate. In this case, the universality class for static properties is fully determined by two exponents. 
Dynamics must also be considered so that one introduces a dynamical critical exponent. 

We consider the case where the transition is at the antiferromagnetic wave vector in d dimensions: Q2= (tt, n) , 
Q3= (it, 7T, it) etc. Since Qd is at the corner of the Brillouin zone, the spin susceptibility \o (Qd) is always, by 
symmetry, an extremum. This extremum is the absolute maximum at half-filling not only in the nearest-neighbor 
hopping model, but also in more general models with next-nearest-neighbor hopping || Jl5[ |. The nearest-neighbor 
model is discussed in more details at the end of this section. It has some special features resulting from the additional 
nesting symmetry. In the two-dimensional case, we also comment on peculiarities of nesting and on quantum-critical 
behavior |] |]. 



A. Renormalized classical regime. 



As one decreases the temperature sufficiently close to the phase transition, there appears a small energy scale 6U 
that measures the proximity to the phase transition as determined by the Stoner criterion. This scale is defined more 
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precisely in Eq.(|T^). The key physical point is that this energy scale is the smallest. In particular, it is smaller than 
the temperature 

SU<T (9) 

so that the zero Matsubara frequency representing classical behavior dominates all others. The self-consistency 
conditions Eqs.(Q)(^) then lead to a strong temperature dependence of 5U. This is the renormalized-classical regime 
jij. In this regime, the antiferromagnetic correlation length £ becomes so large that |l3| 

£ » fa (io) 

where 

««« = —f ( U > 

is the single-particle thermal de Broglie wavelength and (vp) is the Fermi velocity averaged over the Fermi surface. 
This provides a partial justification for the usual procedure [9 that eliminates completely the Fermionic variables 
and describes the system in terms of collective Bosonic variables, as is usually done in Hubbard-Stratonovich types 
of approaches. 0] pj 

We first show that when most of the temperature dependence of the susceptibility comes from the temperature 
dependence of 5U, the RPA-like form that we have implies that in any dimension the dynamical exponent is z = 2 
while the classical exponent j/i> = 2 — n takes the value j/u = 2. The other classical exponent v is determined from 
the self-consistency condition Eq.(Q). We show that the corresponding universality class is the same as the n — > oo 
limit of the O (n) classical model. This universality class is known in turn to be the same as that of the spherical 
model. jl6| We conclude this discussion with the lower critical dimension d = 2. There the exponent v cannot strictly 
be defined since, as was shown before |7j], the correlation length diverges exponentially at zero temperature instead 
of diverging as a power law at finite temperature. This behavior is also the one expected from the n — > oo model, 
although nesting leads to different temperature dependences that are explained further. 



1. Exponents j/v and z in arbitrary dimension 



The antiferromagnetic transition is characterized by the appearance of a small energy scale, or equivalently a large 
correlation length, in the retarded spin susceptibility 



xf„(q,^) 



i-i[/ spX £(q,o;)' 



The small energy scale is set by 

SU = U m f,c - u sp 

where the temperature-dependent "mean-field critical" interaction 

C/ m/ , c = 2/ X o(Qd,0) 



(12) 



(13) 



(14) 



is the temperature-dependent value of U sp at which a phase transition would occur according to mean-field theory. In 
the vicinity of this point the small energy scale SU allows us to approximate x2>(q> w ) by expanding the denominator 
near q « Q d and w~0to obtain, 



x™(q + Q d .w) «C 5 



1 



i + q 2 £ 2 ~iut 2 /D 



where the antiferromagnetic correlation length is defined by 



w 



1/2 



(15) 



(16) 



with the microscopic length scale set by 
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*o 2 



d 2 Xo (q,0) 



2X0 (Qd) dql 



The microscopic diffusion constant D is defined on the other hand by 



1 

D 



TO 
?0 



where the microscopic relaxation time is, 



To 



f 



xo (Qd) 



(17) 



(18) 



(19) 



This relaxation-time is non-zero in both models B and C where the Fermi surface intersects the magnetic Brillouin 
zone. 

In the presence of a large correlation length £ the scaling q ~ £ _1 and a; ~ £~ 2 justifies the neglect of higher-order 
terms in the expansion Eq.([T5|). Comparing the approximate form Eq.(|i"5|) with the general scaling expression 



where X (q£, lu^ z ) is a scaling function, we immediately have the announced results, 

z = 2. 



1 = 2 



(20) 



(21) 



The Fisher scaling law n = 2 — ^ shows that the anomalous exponent r\ vanishes as in mean-field theory. In the 
following paragraphs, we compute the remaining exponent v to show that above four dimensions we do recover 
mean-field theory v = 1/2 while for 2 < d < 4, we have the n — > 00 result v = 1/ (d — 2). 



,2. Exponent v in 2 < d < 4 and equivalence to spherical (n — > 00) model. 

The correlation length exponent is determined by solving self-consistently Eqs.(Q) and (||) for the quantity (n^ni) = 
Usp/U . In general, we do this numerically using some technical tricks that are discussed below. With this procedure, 
no arbitrariness is left in the cutoffs, that are entirely determined from the microscopic Hubbard model. However, to 
study analytically the critical behavior, we notice that there is a crossover temperature Tx below which the presence 
of the small energy scale SU << T makes the zero Matsubara frequency component in the sum rule Eq.(|7J) much 
larger than all the others. This is the renormalized classical regime discussed above. Its existence is a manifestation 
of critical-slowing down, u> ~ ~ SU near a phase transition. Using the approximate Lorentzian form Eq.(|l5|) for 
the iu> = iq n = component we rewrite Eq.(j^) as follows, after a trivial shift of integration variables, 

2 _ 2T f d d q e (99 , 

In this equation, 

ct 2 = n - 2(n T n x ) -C < 1 (23) 

is the local moment n — 2(n^n{) minus corrections C that come from the sum over non-zero Matsubara frequencies 
(quantum effects) and from the terms neglected in the Lorentzian approximation, namely those coming from short 
distances (q — Q) 3> £~ 2 . 

Contrary to the strong-coupling case, and contrary to more usual approaches 0] , a 1 here is temperature dependent 
because both (n-j-n^) and C are. Nevertheless, to find the critical behavior analytically, it suffices to notice that 
this dependence is regular. In fact, we have that when T ~ Tx, the double occupancy can be approximated by 
(n^ni) = U S p/U w U m f :C /U when SU — > 0. 

At the Neel temperature, Tjv, the correlation length diverges, £ = 00, so that Eq.(p2|) determines the Neel temper- 
ature through 



„ 2 2T N f d d q 1 

a = 



U sp e J (2tt)<V 



(24) 
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The wave vector integration is cutoff at large q by the Brillouin zone (— tt < qi < n for any component Qi) so that 
the only divergence occurs from g = 0ind<2(g = 0— > Q = Qd in the original integration variables). Since the 
left-hand side of this Eq.(^4|) is finite, this divergence prevents the existence of a finite-temperature antiferromagnetic 
phase transition in two dimensions or less. 

To find the correlation length exponent in 2 < d < 4, one rewrites Eq.(B3) in the form, 



a 2 = 



2T 



d d q 



U sp e J (2ir)d 



e 



1 + q 2 ^ 



r 



2T 



d d q 



Using the expression for the Neel temperature Eq.(E4|), this last expression becomes, 



a 2 



1-^1 = 

Tn 



Since the integral converges at q£ 
one finds, 



2T 

TJ W 



t 2-d 



d d (gfl 
(27r) d 



1 



1 



1 + q 2 ? q 2 £, 2 



(25) 



(26) 



oo for 2 < d < 4, it can be replaced by a ^-independent negative number and 



1 ' - 1 



l/(2-d) 



T 



- 1 



(27) 



which gives, 



1 



d-2 



(28) 



This exponent and 7 jv = 2, found in the previous section, are the same as the one for the Berlin-Kac spherical model 
ijilif or equivalently for the generalized Heisenberg model where spins are n components vectors and n — > 00. In three 
dimensions, this leads to 
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77 = and (5 = 5. 



(29) 



For comparison, numerical results |l7[ for the 3D Heisenberg (n = 3) model give v ~ 0.7 and 7 ~ 1.4. 

Above d — A, one recovers the mean-field results 7 = 2z; and v = 1/2. This last result follows from the fact that in 

d > 4, the integral in Eq.([2^) is dominated by the large momentum cutoff so that for £ >> 1, I 1 — j ~ £~ 2 J d d q/q i . 



3. Two-dimensional case 



We have already proven in the last subsection that the transition temperature vanishes in two dimensions. The 
correlation length may be found in the renormalized classical regime directly by performing the integral Eq.J^) 
in d = 2, 

£ = Co (U sp /5U)^ ~ A^ 1 exp(7ra 2 ^[/ sp /T) (30) 

where A ~ -k is usually of the order of the size of the Brillouin zone, but not always as we discuss below. 

In d — 2, we call Tx the temperature at which SU is much smaller than temperature and the magnetic correlation 
length £ grows exponentially. While in higher dimensions a phase transition occurs at finite temperature, in d = 2 
the critical regime with an exponentially increasing £ extends all the way to zero temperature. For example, the 
temperature Tx isplotted as a function of filling in the two-dimensional nearest-neighbor Hubbard model for U = 2.5 
in Fig.l of Ref. [Q. In this reference, Tx is called a quasi-critical temperature. We stress that there is a range 
of fillings near half-filling where at Tx it is the antiferromagnetic wave vector that grows, despite the fact that at 
zero-temperature the phase transition would be at an incommensurate wave vector. 

The exponential growth of the two-dimensional £ clearly suggests that small 3D effects existing in real systems 
may stabilize long-range order at Qd=3, before T = 0. We later characterize the crossover driven by a small 3D 
hopping parameter tu <C tj_ from two-dimensional critical behavior to three-dimensional critical behavior. But 
before, we comment on the two-dimensional quantum-critical regime and on peculiarities induced by nesting in the 
renormalized-classical regime. 
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B. Quantum-critical regime 



When there is a critical value of the interaction U c at zero temperature where one finds a paramagnet for U < U c 
and an antifcrromagnet for U > U c , then the T = 0, U = U c point of the phase diagram is a quantum critical point. 
The vicinity of this point in two dimensions has been studied again recently [2|. In order to study such a regime 
within the Hubbard model at half-filling, one must introduce next-nearest-neighbor hopping since U c (T = 0) = 
at this filling in the nearest-neighbor model. One finds that the TPSC approach has precisely the n — > oo model 
A or model B quantum critical behavior [Q], depending on the specific microscopic model. In particular, £ scales 
as 1/T as one approaches the two-dimensional quantum critical point from finite temperature. Again, in the TPSC 
approach the cutoffs are specified without ambiguity. Model C, the perfect nesting case, is relevant only to the 
renormalized-classical case, as we now discuss. 



C. Peculiarities induced by perfect nesting in the renormalized-classical regime, especially in d — 2. 

The dispersion relation of the nearest-neighbor Hubbard model on hypercubic lattices in arbitrary dimension satisfies 
£k+Q d = — £k- Furthermore, at half-filling the particle-hole symmetry implies that the Fermi surface is fully nested, 
namely \x = so that the equality ek+Q d — A* = — (ek — is satisfied for all wave vectors k. Slightly away from 
half-filling, nesting in the form £k+Q d — t 1 ~ — (ek — A*) is a l so a good approximation at finite temperature as long as 
T > /i, as discussed above. When there is perfect nesting, the zero-temperature critical interaction vanishes (U c = 0). 
Hence the fully nested Fermi surface, referred to as Model C above, does not have the simple quantum-critical point 
described in the previous sub-section. 

When there is perfect nesting, the microscopic interaction-independent quantities £q an d tq have a peculiar tem- 
perature dependence. This occurs because they are derivatives of the susceptibility which itself contains logarithmic 
singularities in the zero-temperature limit. These quantities are evaluated in two dimensions and in the quasi two- 
dimensional case in Appendix A. Dimensional arguments that follow simply from this appendix show that in d > 2 

^~l/(T a lnT- 1 ) (31) 



ro-l/tTlnT- 1 ). (32) 

In d = 2, the lnT -1 is replaced by In 2 T -1 . g§ 

By contrast, in the case of second-neighbor hopping, nesting is lost and the above quantities are temperature 
independent for a wide range of values of the second- neighbor hopping constant. The above temperature dependencies 
are then a special property of nesting. In d > 2 however, the above temperature dependencies are completely negligible 
in the critical regime since near the phase transition one can replace T in the above expressions by Tn- 

The only issue then is in two dimensions where the phase transition occurs at zero temperature. Even neglecting 
logarithms for the moment, one sees that since £ 2 scales as 1/T 2 over a wide temperature range the correlation length 
in Eq.(|30|) scales as exp(csi/T 3 ). By contrast, in strong coupling, or in the non-nesting case of the weak-coupling 
limit, the correlation length scales as exp(cst/T). 

The exp (est/ (T 3 ln 2 T -1 )) behavior is however largely an unsolved problem. Indeed, in the critical regime in two 
dimensions, fluctuations remove the quasiparticle peak and replace it by precursors of the antiferromagnetic bands, 
as shown in Ref. It is possible then that, in this regime, a more self-consistent treatment would lead to £ 2 

independent of temperature, as in the strong coupling case or the non-nested weak-coupling case. It is also likely that 
there will be an intermediate temperature range where the exp (est/ (TWT- 1 )) regime prevails, even if deep in the 
critical regime self-consistency leads to exp(csi/T) behavior. 

It is important to recall that in practical calculations in the TPSC approach, one obtains a numerical value for the 
correlation length without adjustable parameter. For example in Fig.l we present the temperature dependence of the 
correlation length for the two-dimensional nearest-neighbor Hubbard model. As discussed in Appendix A, in this case 

Co 2 ~ 0.021U mf , c tiai/T 2 (33) 

and U sp ~ C/ m /,c so that from the slope of the plot and from Eq.(|30|) one finds a 2 ~ 0.21. From the plot we can also 
extract A -1 ~ 0.022 so that £ is known without adjustable parameter. Appendix B explains physically the orders 
of magnitude taken by a 1 and A" 1 in this model. Similar calculations can be done for arbitrary band structure. In 
strong-coupling calculations, |L) || one obtains £ ~ A -1 exp(27rps/T) with ps a cutoff-dependent quantity that can 
be evaluated only with Monte Carlo simulations. 
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Another consequence of the temperature behavior of £o m Eq. ( pi| ) is that above Tx there is a range of temperatures 
for which the antiferromagnetic correlation length scales as £ ~ £o ~ 1/T. This behavior should not be confused 
with quantum-critical behavior, even though the power-law scaling of the correlation length is the same. Indeed, one 
finds that the argument of the exponential in Eq.(^o|) is larger than unity in the corresponding regime while in the 
quantum-critical regime the argument of the exponential should be much less than unity. In fact the temperature 
dependence of the staggered susceptibility for T > Tx is also-different from the quantum critical result. 



IV. QUASI TWO-DIMENSIONAL SYSTEMS: RENORMALIZED CLASSICAL CROSSOVER FROM D = 2 

TO D = 3. 

The general discussion of universality in the renormalized-classical crossover from d = 2 to d = 3 appears in 
Appendices C and D. In the present section, we first clarify the various regimes of crossover, according to whether 
or not single-particle coherence in the third dimension is established before the phase transition. Then, we go on 
to discuss the case fii <C T/v < Tx where the SDW instability occurs before interplane single-particle coherence is 
established. More specifically, we find the scaling of the Neel temperature with t±/t\\ as well as the size of the three- 
dimensional critical regime with the corresponding exponents, showing that the results are those of the n — > oo limit. 
We restore the lattice spacing units an along the three-dimensional axis and a±_ in the planes. We assume however 
that the ratio a\\/a± is usually of order unity and numerical calculations are done for an /ax = 1 



A. One-particle and two-particle crossover from d = 2 to d = 3. 

We consider in this section the highly anisotropic situation where hopping between planes tu is much smaller than 
in-plane hopping, t±, 

t\\ <C t± (34) 

as might occur in the high-temperature superconductor parent compound LaiCuO^. In this case, we have that the 
three-dimensional transition temperature to long-range order Tn is always less than the crossover temperature Tx to 
the characteristic exponential behavior of the correlation length in two dimensions 

T N <T X . (35) 

This is so because the microscopic in-plane and out of plane £g lengths satisfy £q- £q . 

The crossover temperature to two-dimensional behavior for itinerant antiferromagnet always satisfies 

T x <t±. (36) 

Two limiting cases are then possible, depending on interaction and on hopping parallel to the three-dimensional axis 

a) Weak coupling or small anisotropy limit: 

T N <T x <$:t [l . (37) 

In this case, when the transition to three-dimensional behavior occurs the three-dimensional Fermi surface is relevant 
since the thermal de Broglie wavelength Vp/(nT) ~ tiian/T is larger than the distance between planes. In other 
words, the three-dimensional band structure is relevant to the behavior of single-particle propagators in Matsubara 
frequencies before the phase transition occurs. Fermions are quantum-mechanically coherent over more than a single 
plane and nesting generally plays a role in the value of the ordering wave vector. The crossover from two to three 
dimensional critical behavior would occur in a manner analogous to the anisotropic Heisenberg model [fl9f |20| ]2l| | 

b) Intermediate coupling or very large anisotropy: 

t\\<.T N <Tx. (38) 

Here, long-range order is established before the single-particle coherence occurs between planes. A phase transition 
occurs only because of two-particle (or particle-hole) coherent hopping. When the phase transition occurs, thermal 
fluctuations are still large enough that coherent single-particle band motion in the parallel direction has not occured 
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yet. There are several ways of explaining physically what this last statement means. For example, it is clear that 
when tii <C T, features of the band structure in the parallel direction are irrelevant for single-particle properties since 
the first Matsubara frequency is larger than the bandwidth in that direction. The motion between planes is still 
in that sense quasi-classical when the phase transition occurs. Another way of saying what this means is that the 
thermally induced uncertainty in the parallel wave vector is equal to the extent of the Brillouin zone in that direction, 
corresponding, via the uncertainty principle, to a confinement within each plane. 

We do not discuss the intermediate case T/v < t\\ < Tx but concentrate instead on the very large anisotropy - 
intermediate coupling limit just introduced. 



B. Numerical solutions 



The scaling behavior of the Neel temperature T/v and of the three-dimensional crossover temperature are derived 
in the following two sections. We first present the numerical results obtained from the solution of the self-consistency 
relations Eqs.(0) and (^). The numerical integration in Eq.(Q) is made possible by rewriting this equation in the form, 

n - 2(n ]ni ) = Toyoi / -£|r ]T [ Xsp (q,i<Zn) - xTp MR.o] + Ta\\a a x f -^xTp M) . (39) 
J (2tt) iqn J (2tt) 

The sum over large Matsubara frequencies can be approximated by an integral in a controlled manner. The subtraction 
in the first integral removes singularities of the integrand and makes the integral well behaved. Since the transition 
occurs in the single-particle incoherent regime Eq.(pq), the integrand in square brackets is independent of q» above 
the Neel temperature. All quantities involving in come from the second integral over the asymptotic expression for 
the susceptibility 

X-(q+Q3,0)^^ 1 + ^ +gM (40) 



where 



^^ti(U sp /6Uy/ 2 , (41) 



£n^(iV^) 1/2 - 



(42) 



To have sufficient precision for large two-dimensional correlation lengths, it is important to evaluate analytically 
£q, as well as the integral of Eq. (|40|) appearing in the consistency equation Eq.(|39"|). This is done respectively in 

Appendices A and C. To perform the second derivatives in the definition of £q~, £q, we expand Xo(Qs) m pow ers of 
t||/T, keeping only the first non-zero term: thus does not differ from the one already presented in Eq. (|33|). It is 
shown in Appendix A that over a wide range of temperatures we have 



and 



d 2 Xo (Q 3 ) 
dq\ 



f2 
T 2 



0{tJT) 



(43) 



d 2 X o (Qa) 
dq 2 



t 

J! 2 



^ + 0(tl/T 3 ) 



(44) 



The inter-plane hopping t\\ in Eqs.(p9|) and (|40D occurs explicitly only in and the above results imply that, 



1 



tj_aj_ 



(45) 



We present numerical results for the nearest-neighbor Hubbard model in units where an = aj_ = 1 and t± = 1. 
The value of TV (t\\) appears in Figs. 2 and 3 for U — 4. In Fig. 2, we clearly see that Tjy becomes almost equal to 
Tx ~ 0.2 for f|| still quite small. The scaling of TV (t\\) shown in Fig. 3 is explained in the following subsection. Fig. 4 



9 



shows the variation of the in-plane correlation length £,± as a function of temperature for various tii, again for U = 4. 
For the purely two-dimensional case tu = 0, one can observe for T < Tx — 0.2 the exponential behavior mentioned 
in the preceding section. At a temperature about 0.16, the in-plane correlation length £j_ is already as much as 10 3 
(in units where lattice space is unity). For 10~ 3 < fii < 10 _1 , £± diverges at the Neel temperature located in the 
narrow range 0.16 < TV < 0.2. For t\\ ^ 0, it is also clear that the crossover to three-dimensional behavior occurs in 
an extremely narrow temperature range. This is explained below. Note that the last few curves on the right-hand 
side are at the limit of validity of our approximations. 

We note that for U — At, we find that at T — TV the local moment is equal to three quarters of the full moment in 
the atomic limit, i.e. n — 2(rtj)(nj)<?j|(0) = 0.75, <?jj,(0) = (n-\n^) / (n^)(ni). This number is only weakly dependent 
on temperature in the range studied. 



C. Dependence of the Neel temperature Tat on tn. Crossover exponent. 



From the discussion of the previous section, Eqs.(p9[)(|4C|)(|22|) we know that the singular part of the self-consistency 
condition may be written in the following form in the quasi-two dimensional case, 



a 2 



2Ta||a 2 L 



d 3 



q 



i 



(27T) 



it)' 



(46) 



The integral can be done exactly, as in Appendix C, and all the results obtained in this subsection and the following 
one can be obtained from limiting cases of this general analytical result, as shown in Appendix D. Here we make 
approximations directly on the integrals since this makes the Physics of the results more transparent. Although 
arbitrary cutoffs appear in the analytical expressions, we re-emphasize that in the numerical calculations of the 
previous section the cutoffs are simply given by the Brillouin zone and there is no arbitrary scale in the results. 

At TV, we have £J 2 = 0. Furthermore, from Eqs.([4i"|) and ( [i"2| ) we find = £o/£o~ so that tlie- above integral 

Eq.([46"|) takes the form 



a 2 



2T N a\\a 2 L f d 2 q± 



dq\\ 



2?r 



(47) 



Using the mean-value theorem for the integral over qu we have 



a 2 



2T N a\ 



fq± 



t ii 



A 2 



(48) 



where A is a constant that we need not specify. It is contained in the range < 



the same as the one that determines the correlation length in two dimensions Eq.(22), hence at Tn we have that 



A 



< 7r/ai|. The above integral is 



ei 2 {Tj 



N) 




A 2 oc 



*ll a ll 

1 1 a i 



A 2 . 



(49) 



Comparing with the general theory of Appendix D where it is argued that £x (Tjv) ~ ~k\ , we see that <f>/v = 2. In 
other words, <j)/v — j/u — 2 and the crossover exponent |l^] is here equal to 7 as is usually the case in the n — > 00 
model. p3| We obtain, using the expression (|30|) for the correlation length in two dimensions, 



1 



Inf- 
ill 



(50) 



where c is a non-universal constant of order unity. 
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In the special case of perfect nesting (half-filled nearest-neighbor hopping model), the microscopic length £q- is 
temperature dependent, as shown in Appendix A 

,,,,2 2 0.085 1 

fe) ~ ai — ^y- (51) 

Using this result as well as f/ m /, c = 2/x (Q2) ~ £4p at ^Tx ~ ^V gives the scaling illustrated in Fig. 3 for the case 
ox/011 = 1 namely 



1 T 2 

1 J AT 



m/,c 



In J- 



(52) 



The logarithmic behavior in Eq.(|50|) is typical of systems that undergo a dimensional crossover from their lower 
critical dimension. For example, the analog of Eq.(^9|) in the anisotropic Heisenberg case would read, JlS| 

G a (Tn) ~ (i) A 2 (53) 

leading to [[l9| T^ 1 ~ In Cjj^j • The above results Eqs.(^9|) and ( ^3|) are suggested by the simple RPA-like form 
X3d ~ X2d/ (l - J\\X2d) with J\\X2d ~ 1 at the transition and \2d ~ ~ (]_ ~ ex P (Jxcst/T). As in the previous 
section, the quantity [7 sp f^o"'"j pl avs a r °l e analogous to the exchange constants J±\\ r In the perfect nesting case, 

these effective exchange constants would be temperature dependent since U sp fCo*'") ~ U S p (t±,\\) 2 /T 2 . Note that 
in the crossover from one-dimensional Luttinger liquid behavior to three-dimensional long-range order, the effective 
exchange constant J± also scales |h| |ll]] as ut\jT 2 , with u a running coupling constant. The one-dimensional Fermi 
surface is always nested. 



D. Size of the three-dimensional critical region. 

The singular temperature dependence of the correlation length is obtained from the equation 

^ 2 2Ta\\a 2 L f d 3 q 1 



SU J (2nyi+tiql+qql ^ 

Since the ratio = Co/Ccf i s temperature independent, a simple change of integration variables shows that near 

Tat the scaling of both correlation lengths with temperature is identical to the isotropic three-dimensional case. In 
other words, the critical behavior near the phase transition is that of the three-dimensional system. However, as one 
increases the temperature away from T/v, the correlation lengths can decrease until £11 <C an while at the same time 
3> a±. When £|| <C an, the integral Eq.(|5^) is essentially two-dimensional and for ^> a± one should observe the 
characteristic exponential temperature dependence of the two-dimensional correlation length. 
As usual the definition of crossover contains some arbitrariness, so let us choose 

£|| Cn = oil (55) 

as the definition of the crossover temperature T* between d = 2 and d = 3 critical behavior. In that regime, the 
correlation length Eq.(|55|) scales with temperature as in the d = 2 regime Eq. (^0|) except that, as argued before, £|| is 

smaller by a factor Ko/£(j~^ = (^ll a ll) / (^-i- a -i-)i hence we obtain for T* 



1 



7T0" 



'U (f 1 



. 2 



(56) 



where c' is a non-universal constant of order unity. The size of the crossover region is thus 



1-^ = 7V al 2 {c-J)= {c r^ ^r- (57) 
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The above results Eqs.(|55|) to (57) are as expected from the usual theory of critical phenomena exposed in Appendix 
D. In particular, the scaling of T* with t\\/t± is the same as that of T/v. The smallness of the crossover region from 
d = 2 to d = 3 critical behavior in Fig. 4 follows from the above considerations. The smaller is T/v, the smaller is T*. 
The above situation should be contrasted with the problem of crossover from d = 3 to d = 2 in Helium films, studied 
by Fisher and Barber Bjj ] many years ago. In that case, power law scaling occurred everywhere, giving quite different 
expressions for the scaling of T* and T/v . 

Given (T*) = fjj (T*) {U sp /6U (T*)) 1/2 , and £ ' - t h the above relation £,| (T*) = a\\ means that SU (T*) should 
scale as t|. Similarly we should have SU (T/v) ~ ijj. We checked numerically [^5| that the scaling with iy holds for 
t|| < 0.05 in the half-filled nearest-neighbor model with U = 4tj_. 



V. CONCLUSION 



We have shown that the TPSC approach allows one to study all aspects of nearly antiferromagnetic itinerant 
electrons in one-band Hubbard models. The method is in quantitative agreement with Monte Carlo simulations in the 
non-critical regime M Q while in the critical regime, (renormalized classical or quantum critical) the relatively weak 
temperature dependence of the local moment leads to the same critical behavior as strong-coupling models to leading 
order in the 1/n expansion, namely in the n —> oo limit. There is no arbitrary cutoff so that all results can be obtained 
as a function of lattice spacing, hopping integral and interaction parameter. Fermi surface effects are apparent, in 
particular in the case of perfect nesting where the two-dimensional renormalized-classical correlation length diverges 
as exp (est/ (T 3 In 2 T)) instead of exp (cst/T). 

We have applied the method to a detailed study of the renormalized-classical crossover from two to three dimen- 
sions where we have highlighted the existence of a regime where the three-dimensional Neel instability occurs before 
thermal fluctuations become small enough to allow coherent sing le-p article band motion between planes. An anal- 
ogous phenomenon occurs in quasi-one dimensional systems fic| |llf| . The TPSC approach can be applied to study 
realistic cases. For LaiCuO\ we will show in a subsequent publication that with second-neighbor hopping one can fit 
experiments on the magnetic structure factor. 

The generalization of the TPSC approach beyond leading order in 1/n is left open. Also, the effect of self-energy 
feedback Jl3| on exp (est/ (T 3 ln 2 T)) behavior of the correlation length in the two-dimensional nesting case should 
be cleared in further studies. Finally, the universal d = 2 to d = 3 crossover discussed in Appendix D should be 
investigated beyond leading order in 

We are indebted to C. Bourbonnais for numerous discussions and key ideas. We also thank David Senechal for 
discussions and D.S. Fisher for pointing out Ref. p^|. We acknowledge the support of the Natural Sciences and 
Engineering Research Council of Canada (NSERC), the Fonds pour la formation de chercheurs et l'aide a la recherche 
from the Government of Quebec (FCAR) and (A.-M.S.T.) the Canadian Institute for Advanced Research (CIAR) 
and the Killam foundation. 



APPENDIX A: Q ' x AND to IN THE CASE OF NESTING 



In this appendix we derive expressions for the out of plane £ ' and in-plane £j" microscopic lengths 



11,-L _ 



-1 



d 2 X0 (q,0) 



2xo (Qd) dqf u 

as well as for the microscopic relaxation time To in Eq.(|l9|) 

1 d X $ (Q d ,u>) 



TO = 



Xa (Qd 



diuj 



UJ = 



(Al) 



(A2) 



for the quasi two-dimensional antiferromagnet Q3 = (tt,tt, it) , in the regime ty <C Tn < Tx of Eq. (|3q) . We also 
assume that \x -C T so that the maximum of the static susceptibility is at Q3 even away from half-filling. 
We start from the retarded Lindhard function in d dimensions 



d d k /(£k+q ~ M) ~ /(ek - M) 
l;z ( 27T ) d u + in- e k+q + £k 



(A3) 
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where / is the Fermi function, /i the chemical potential (fi = at half-filling for our Hamiltonian but the expressions 
quoted here are more general than for the half- filled case). For nearest-neighbor hopping, we have the nesting property 



e k +Q d = -£k (A4) 

that can be used to rewrite 

Xo (Q* = 2 J dEN d (E) u + ir} + 2E (A5) 

where Nd(E) is the single-spin density of states for the given dimension. 
In the limit /i « T we have for the static susceptibility 

xo(Qd) = xf (Qs, o) = 2 1 d£jv d (£) 1 ~l f E [E) (A6) 

so that in two dimensions Xo(Q2) ~ In 2 (i/T) while in three dimensions Xo(Qa) ~ ln(t/T). In the quasi two- 
dimensional case with in <§C T the two-dimensional value of xo is an accurate approximation. The numerical values 
of Xo(Qd) are in practice easy to obtain from numerical integrations. 
For the microscopic relaxation time when /i « T we start from 

Im X *(Q rf ,u;) = 7riV d (|)tanh (^) . (A7) 

In two dimensions, the logarithmic divergence of the density of states Nd(^) at the van Hove singularity makes the 
zero-frequency limit of the microscopic relaxation time ill-defined. Nevertheless, van Hove singularities are usually 
washed out by lifetime effects in more self-consistent treatments so that one expects that for lu << T one has 
dxo (Qd,w) /du)\^_ ~ 1/T leading to the temperature scaling of To ~ (TlnT -1 ) in d > 2 described in the text. 

We move on to evaluate analytically the wave vector derivatives in the regime t|| <C TV < Tx- Keeping for a while 
a general notation where i is some direction (x,y, or z), one can write 

d 2 Xo n o 2 f d 3 k d 2 C . s t f d 3 k dC 

-^- = -8t 2 a 2 . , 3a2 sm 2 {h + q t - 4Ua t / -—3- cos(fci + ft), A8 

»9i Jbz (27r) ci <9e£ + 7 SZ (27r) d <9e k+q 



where 



L-lEk+qjekJ — 



£k+q — £k 

Assuming tu/T <C 1, we evaluate second derivatives to the lowest non-zero term in powers of t«/T. For = q», the 
leading term in Eq. (A8) gives a t| contribution if we keep t\\ = in the integrand. The second term gives also to 



leading order a quadratic contribution in t\\. The spread of the Fermi factors over an energy interval of order T allows 
us to neglect all other dependencies in t» and to perform the integral in the third direction trivially, enabling us to 
rewrite the remaining integral in terms of the two-dimensional single-spin density of states N2 (E) . After some algebra 
we get 

92 *°< Q -> —adj f\ ENME){ i^±A±fE^l + w (E } + 0(( | )3) ( A9) 



where /' is the derivative of the Fermi function. Using the expansion of Fermi functions and derivatives near E = 0, 
it can be shown that the integrand in the preceding equation is finite at finite temperature. Indeed as E/T — > 0, it 
behaves as N2D(E)f"'(fj,), where /"' is the third derivative of the Fermi function. At low temperature, approximating 
the integrand by N2d(E) f"' '(/i) over an energy interval T shows immediately that 



d 2 xo(Q 3D ) _ 2 *f 

I J>2 



,2 -«-^7- (AID) 



More precisely, 1/T 2 shou ld b e multiplied by a logarithmic correction that comes from the 2D density of states. 
Numerical integration of (A9) shows that this t 2 -^ behavior occurs on a wider range of temperature than first 



expected: T = 0.2f is already in this regime. 
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get 



We can evaluate the in-plane qi = q± derivative in the same spirit. This time take t\\ = from the start since the 

/ T- VT., i 1,,,. 1 :)2 , (f\ \ /:),.- , r» i /';),,- \ 



leading order is in t\/T 2 > ^/r 2 . We thus have <9 2 Xo(Q3)/<9g 2 ~ <9 2 Xo(Q2)/<9g 2 . After tedious algebra we finally 



d2 ^0^ * jf <*^M£) [\ if'(E + M) + /'(£ - M)] + 



2E 



it 



-/!,] / ( /BM(£)|i[/"(£ + M )+f(£-/ 1 )] 



/'(£ + M )+/'(£- M ) l-/(S + M)-/(^-M) 



£ 2 



where the integral 



M(E) 



dq x A sin q x 



d 2 k 
(2^ 



£ 3 



(All) 



(A12) 



can be interpreted as an average over the surface of constant energy E of the square of the Fermi velocity in the x 
direction times the density of states at this energy. It can be evaluated analytically as 



M(E) = J, { 2E( fc 2 ) - A(l + Jl) F(fc2) + |l n( o 2 , k 2 )}. 



(A13) 



Here F(fc 2 ), E(fc 2 ) and Tl(a 2 ,k 2 ) are complete elliptic integrals of respectively first, second and third kinds, with 
k 2 = 1 - E 2 / [l&t]) , and a 2 = 1 - E/ (U±). Again at E/T = the integrand in Eq. flAlip is well defined, and 
using Fermi function expansion it can be shown that at low temperature d 2 xo(Q)/dq 2 _ scales as a^i^/T 2 , with this 
time no logarithmic prefactor as before. More precisely we found numerically for a wide range of temperature (wider 
than the range studied in the main text) the following behavior 



5 2 X0(Qe 



dq 2 . 



-0.085ai^ 



(A14) 



and correspondingly, 



d 2 Xo (q,0) 



2X0 (Q d ) dqj 



0.085a 



77 j-2 
2 u mf,c 



4 T 2 ' 



From Eqs.(AlO) and (A14) and a numerical evaluation of the corresponding quantities, one finds the following scaling 



h a \\ 

t\a\ 



(A15) 



To conclude this appendix let us stress the fact that expanding Xo( c l ~ Q21)) ^° ^ ne second order using Eq. (All) 
to obtain the asymptotic form of the 2Z?-spin susceptibility is valid as long as the maximum of xo is & t {^i 7r )? which 
is more general than half-filling. Indeed, by symmetry, the first derivative of the free susceptibility at Q2D is zero 
for all fillings n and temperature and, as discussed before [p6| H, at finite temperature and away from half-filling 
the absolute maximum of the free susceptibility can be at (71", it) even if it is not the case at T = 0. This behavior 
can be observed in Figure 5 where the in-plane second derivative is plotted as a function of temperature for various 
values of band filling. When the second derivative goes to zero there is a shift in the wave vector maximizing the free 
susceptibility. Whether the magnetic transition will be commensurate or incommensurate at a given filling depends 
on the interaction U because by changing U one can change the ratio [ijTx and because the nature of the final 
three-dimensional order depends very much at which wave vector correlations start to grow below T x . 

Calculations presented along the above lines do not allow us to study the case where the maximum occurs at an 
incommensurate vector since we need the analytical expressions for the second derivatives of \o to perform very 
accurate numerical calculations (When the wave vector q is different from (tt, tt) we do not have anymore the useful 
simplification: ek+ q = — £k allowing us to replace the (k x , k y ) integration in Eq. ( |A8| ) by a simpler integral on the 2D 
density of states.) Progress is nevertheless possible numerically within the TPSC approach. 
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APPENDIX B: ESTIMATES FOR a 2 AND A 1 IN THE NEARST-NEIGHBOR MODEL. 



In this appendix we provide estimates for a 2 and A -1 in the isotropic two-dimensional nearest-neighbor model. 
The surprisingly low numerical values a 2 ~ 0.21, A" 1 ~ 0.022 obtained for U = 4 in the text are special to the Model 
C perfect nesting case. 

We first rearrange the self-consistency equation Eq.(Q) to isolate the asymptotic behavior, as we did in Eq.(|3^) but 
here in two dimensions and with a± = 1. 

n - 2(n T n x ) = T [ ^xT P M) + T [ ' £ [Xs P (q,^n) - X? P (q,0) S n , ] . (Bl) 

J V 2lT ) J i 2lT ) iq n 

It is usually assumed that the last integral on the right-hand side is weakly temperature dependent and it is included 
with the left-hand side to define a 2 . This procedure usually suffices for reasons we will see below. For a more accurate 
estimate of a 2 close to Tx we use the Euler-Maclaurin formula to approximate the sum over Matsubara frequencies 
larger than the zeroth one by an integral. Recalling also that Xsp (Q^lri) — Xsp (q, —iQn) we have 

n - 2(n m ) = T [ -^xT P M) +T f -^t [ Xsp (q,0) - x% (q,0)] 
J [2tt) J (2ir) 



d 2 q . [°° dX f d 2 q 

(27r) J2ttT 27r J (27r) 



-T / Xs P (q,i?i) + 2 / — / --^Xsp (q,iA) . (B2) 



To recast this result in the same form as the consistency condition Eq.(|22|), we first note that a more satisfactory 
definition of a 2 than the one given in Eq.(23) would be 

r°° d\ f d 2 n 

a 2 = n-2( ni n L )-2 — Xsp (q,iA) . (B3) 



2ttT 



2ir J (2tt) 



Also, the coefficient of the term linear in temperature on the right-hand side of Eq.(22) would not only include the 
asymptotic Lorentzian form but also a correction from the deviation to Lorentzian and another correction from the 
first Matsubara frequency. Overall then, a more accurate expression for the consistency condition is given by the last 
definition of a 2 and 

?2 = T {/ (io 1 ^ (q ' 0) + / (Sp [xsp (q ' 0) ~ xTp (q ' 0)] + / ^f Xsp {q,m) } ■ (B4) 

The rest of this appendix is in two parts. We first estimate the left-hand side of this equation, a 2 , and then we 
estimate the integrals on the right-hand side to obtain A^ 1 . 

To obtain a 2 , one should first notice that at the crossover temperature the local moment n — 2{w\n\) is already 
quite close to its zero-tcmperature value. Taking this as an estimate, we have 

n - 2<n T n L ) = 2 f ^ / -^-^ X s P (q,iA) (B5) 
Jo 2lr J (2n) 



so that, substituting back into Eq.(p33|), we have 



, f 2nT dX f d 2 q 
° = 2 / / Xs P q,*A . (B6) 

To estimate this integral for the case of perfect nesting , we note that singularities of Xsp (q,0) near wave vectors 
q = and q = (tt, it) are integrable singularities. We thus use the mean-value theorem to write, in our dimensionlcss 
units 



/ 



d 2 q 

~2 Xs P (q,«A) ~ Xsp (q typ ,iX) (B7) 



As a representative point, one can take <\t vp = (tt, 0) since it is far from both singularities. Using the trapezoidal rule 
to estimate the frequency integral, one has 
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2ttT 
j ~2t7 



Xs P (qtgj) + Xsp (q tw ,27rT) 



0.19 



(B8) 



whose numerical value follows from results obtained for U = 4, Tx — 0.2, U s , 



mf,c 



Xsp 
Xsp (Qtyp 



,2ttT x ) 



0.60 
0.36. 



(B9) 



The estimated numerical value of a 2 in Eq.(B8) corresponds closely to the value obtained in the text from accurate 
numerical solutions. The fact that <7 2 scales roughly as Tx ~ Tm/,c in very weak coupling, as follows from Eq.(p8[), is 
a significant result since a 2 is also related to the size of the pseudogap between precursors of antiferromagnetic bands, 
as shown in Ref. Jl3| . 

To estimate the value of A -1 , we notice that in the usual consistency condition Eq.(p2|), one keeps only the first 
term on the right-hand side of the more accurate expression Eq.(B4). The effect of the other terms is mimicked by 
using an effective cutoff A that is not equal to the Brillouin zone size, as one might have naively expected. In other 
words, the effective cutoff A may be obtained by requiring that 



A ^X-(q,0) 

Z7T 



(27T)- 



[Xs P (q,0) - xT P (q,0)] 



q 



(2tt) 



2Xs P (q,%) 



(B10) 



When there is no nesting, the quantity £o is relatively small at the crossover temperature, meaning that the asymptotic 
Lorentzian form is not so peaked and should be a good estimate of the susceptibility over much of the Brillouin zone. 
Because of the slow decay of the asymptotic form Xsp (q>0)> the second integral should in fact be negative and should 
partly cancel the last integral so that we should have A ~ jr. By contrast, for perfect nesting £o ~ ^-/T is large, 
as seen in Appendix A, meaning that in this case the asymptotic form Xs P (q>0) is valid only in a narrow range of 
q values. Over most of the Brillouin zone, away from the maximum, the true susceptibility Xsp (q ; 0) is larger than 
the asymptotic one \1 P (q>0) because the latter decays rapidly away from the maximum while the true one has an 
extremum at both the Brillouin zone corner and center. The same arguments as those used to evaluate integrals for 
ct 2 allow us then to estimate 



/d 2 
[Xs P (q,0) - xTp (q,0)] * Xsp (q tOT ,0) a 0.60 

f d 2 q 

/ TZ—2Xsp (q,%) - Xsp (qty P ,2TrT x ) ~ 0.36 
J (27r) 



so that the equation Eq.( BlO ) that determines the cutoff becomes, with U sp ~ 2 and £ 2 (Tx) — 1, 



qdq 2 



U S p£() £ 



2 t-2 



7r exp 



0.96 



0.96 



0.016. 



(Bll) 
(B12) 

(B13) 
(B14) 



Although the difference with the numerically accurate result seems relatively large, one should really compare the 
estimates of In A -1 . The above estimate, ln0.016 ~ —4.1, differs only by roughly 10% from the estimate, lnO.022 = 
—3.8, obtained from a logarithmic plot of the numerically accurate solution. 



APPENDIX C: EXACT RESULT FOR / D 3 Qxsl 



In this appendix we find the integral of the asymptotic part of the spin susceptibility near Q3 
Xas(q + Q3,0) be the approximate spin susceptibility near Q3 obtained in Eq.(f40|) with q\ 



(71", 7T, 71"). Let 

= Qz 



ql + ql and q\ 



First we integrate in the ^-direction from — A|| to A||, with A|| = 7r/aii, then change to polar coordinates in the plane 
and integrate on a circle of radius A^ to finally obtain 



d 3 q 

D (270 



s X^(q + Q 3 ,0) = 



nUap (Co") a || 



A||i| 



-arctgAy^H 
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(CI) 



This analytical result provides another route to obtain the Neel temperature and the d — 2 to d = 3 crossover as 
discussed in the following section. 

APPENDIX D: EXTENDED SCALING HYPOTHESIS AND UNIVERSALITY FOR THE 
RENORMALIZED-CLASSICAL D = 2 TO D = 3 CROSSOVER. 

We first briefly recall the results of Ref. j^jj on universality of crossover scaling functions in anisotropic systems. 
The discussion usually centers on anisotropy in spin space rather than position space but the results are generally 
applicable. Suppose one has a very small anisotropy g. Sufficiently far from the transition, the critical behavior will 
be that of the isotropic fixed point and should be described by the extended scaling hypothesis for the singular part of 
the free energy density. The same extended scaling hypothesis follows for other thermodynamic response functions. 
We use the symbol (w) to mean "asymptotically equal to" and (~) to mean "scales as". Let us concentrate on the 
magnetic susceptibility 

x(g,t)MAL-iX(Bg/t*) (Dl) 

where 

with <p the crossover exponent and T c (0) the value of the transition temperature at zero anisotropy g = 0. It is clearly 



the large value of the correlation length that validates the scaling hypothesis. The scale factors A and B in Eq.(Dl) 
are non-universal, but the scaling function is. The value of A for a given model is fixed by the normalization condition 
*(0) = 1. 

Near the true transition temperature at the anisotropic fixed point, the susceptibility should obey the usual result 

X(g,t)MA(g)i~* (D3) 
where quantities with a dot refer to properties of the anisotropic fixed point, and 

t( r, 8) = fI^W T ° (g >:f° (0) W + M g ). (D4) 



r c (o) J \ t c (o) 



The two expressions for the susceptibility Eqs.( |Dl[ ) and (D3) are consistent only if the crossover scaling function 



X (Bg/t^ is singular as a function of its argument, namely 

lim X (x) = X ( 1 - — (D5) 

where Xq is a universal amplitude while x c is a g and t independent universal number. The definition 

x = Bg/(t(g)f (D6) 
immediately implies that the transition temperature is at 

t e (g) = (Bg/xcf* . (D7) 

The generalization to the d — 2 to d — 3 crossover is not completely trivial because in d = 2 the correlation length 
is not a power law of temperature for O (n) models with n > 1. Fisher and Barber in their study of crossover in 
helium films have considered the case where the system is three-dimensional at high temperature and two-dimensional 
at low temperature, opposite to the situation we consider. Furthermore, the transition temperature is finite in d = 2 
helium films. Kosterlitz and Santos p2| did consider the case of interest here, both within a one-loop renormalization 
group approach, and in the spherical model. To cast the results of the latter paper in the language of the extended 
scaling hypothesis, it suffices to recall the usual hypothesis that the divergence of the correlation length in the plane is 
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at the origin of the scaling behavior. Hence, t (T) v can be replaced everywhere in the above equations by a function 
of absolute temperature that scales with T in the same way as the two-dimensional correlation length |28| 

£ 2d (T) = T Q exp(C*/T). (D8) 

In this expression, we have allowed for a possible algebraic preexponential factor. For example, to one loop order 
p2f in the momentum-shell method [^9) the preexponential factor is a = (n — 2) _1 while to two-loop order |30| as 
well as in the n — ► oo limit, only the exponential is present, a = 0. In addition t o th e non-universal quantities A 



and B defined above, we now have an additional non- universal constant C in Eq.(D8). This is not fundamentally 
different from the usual case where the relation between t and absolute temperature also involves a non-universal 
constant, namely T c (0). The only difference between the itinerant case and the usual n— vector model is that C can 
be temperature dependent in the case of nesting, as discussed in the text and in Appendix A. When there is no nesting 
symmetry, C is temperature independent. In the strong-coupling limit, one usually defines C = 2irps- 
With the above t (T) v — > ^2d (T) hypothesis, the extended scaling hypothesis becomes 



A.sp 



(Q d ,0)^A^ v x[Bg^ v ) (D9) 



where g = (ti\/tx) plays the role of the anisotropy parameter in the case we have considered in detail in the text. 
The function X (x) is a universal function that we normalize to X (0) = 1. With precisely the same asymptotic form 
as in Eq. ( p5| ) , simple power series expansion in powers of T — Tn allows one to recover the correct critical behavior 
near the three-dimensional Neel temperature. Hen ce, the Neel temperature is given by x c = B (tn/t±) t^d" so that 



with the n — ► oo result <p/v = j/i> = 2 and Eq.(D8) one recovers the result of the main text 

with C taking its appropriate temperature-dependent value in the nesting case. 

We conclude by an explicit calculation of the universal crossover function for the staggered susceptibility in the 
n — » oo limit. In this case, a = in Eq.(|D8|). The general form of the susceptibility is given by Eq.([l5|) with 

e - a (u sp /5u) 

XS(Q*0) w Jj- (Dll) 
The value of SU is in turn obtained by solving the self-consistency condition Eq.(|3^ 

a 2 =T ail alj J^ x -(q,0) (D12) 
where a 2 takes essentially its d = 2 value with very small corrections. We can use the result of the previous appendix 



Eq.(Cl) for the integral. It is then convenient to rewrite the result of the integral in terms of the following dimensionless 



variables 



a = TT7T = . V ,2 = cst T 1 ( D1 3) 



Aiei Ai(^) 2 " v 



u = A„eu =A\\Q(U ap /SU) 



-1/2 



(D14) 



Since we assume that we are in the scaling regime, namely the one where the two-dimensional correlation length is 
very large, (Aj_£j_) » 1, we can use 



OL XL 



to expand Eq.(Cl) and write 



lna" 1 / 2 +ln 



u 2 + 1 



1/2 



1 arctan (u) . 

u 



(D15) 



(D16) 
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If o ne so lves the above implicit equation for u, then the susceptibility Eq.(Dll) follows immediately from u 2 in 
Eq.( |D14| ) since 



Af(dj U sp 

Note that u is a function of the dimensionless quantity x defined by 



2u 2 

x'J P (Qdfi) = — . (D17) 



x = a£% d = a exp (2C/T) = ak\^\ (D18) 
C = — p y ° 1 (D19) 



as may be seen by exponentiating the implicit equation Eq.(DlE) 

,2 X ( 2 A 

exp 2 arctanu I . (D20) 

\ u ) 

Before explicitly solving this equation in limiting cases, let us express the universal scaling function X in terms of 
u. The last equation for the susceptibility Eq.( D17j ) may be rewritten with the above definitions as 

X£(Q*0) = 7^M = A&X (x) = - M X (x) (D21) 



2 . 

sp 



where 



X ( x ) = (D22) 



and 



A= 5 . (D23) 

Ai (C X ) U sp 

The universal crossover function X (x) is plotted in Fig. 6. Let us check various limiting forms analytically. This 
will allow us to recover all the cases studied in the main body of the paper. First, the two-dimensional limit is the 
one where u — > 0. In this limit, the implicit equation ( p2C| ) reduces to x = u 2 . This verifies that we have the proper 
normalization X (x) = 1. The three-dimensional limit is the limit where u 2 — > oo. In this limit, 

x c = e 2 . (D24) 

Keeping the next term in the 1/u expansion, we have 

lim X (x) = (-Y ( 1 - (D25) 

x— >x c \eJ \ x c J 

hence the universal constant X takes the value (7r/e) 2 . As expected the susceptibility exponent in d = 3 is 7 = 2. 
The Neel temperature follows from 

x c = aexp(2C/T A r) (D26) 



T N n >v Audi 



(D27) 



Finally, the d — 2 to d — 3 crossover temperature is given by u = Aii^ii = 1. Obviously, the criterion A||£|| = 1 is 
subjective. We could take A|| ^|| to be equal to any other finite number. For definiteness however, we continue with 
A||£|i = 1. Substituting in Eq.( D20 ) we have x* = \ exp (2 — hence the crossover temperature T* is given by 
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aexp(2C/T*) 



(D28) 



Comparing with the equation for the Neel temperature Eq.( D26 ) we find that the scaling of exp(C/T*) with the 
anisotrop y p arameter a is the same as that of exp(C/Tjv). More specifically, we find, in agreement with the main 
text, Eq.(j57|), that the size of the crossover region is given by, 



Tn 



Tn In (x c /x*) = \ = log C2d(Tjv) (x c /x*) . 



2C 



N) 



(D29) 



In other words, the size of the crossover region, calculated in red uce d temperature, decreases with Tn- 
To complete the relation with the general functional form Eq.flM) postulated above, note that if 



ti 



(D30) 



then x — Bg£$l" — ott^d with <p/v = 2 implies that B = a/g is a number 



B 



A IMI 



(D31) 



that is independent of g because of the scaling A|£|/ ^ 2 ±£,± ~ (t\\/t±) that follows from Appendix A, Eq.( A15 ). 

The universal crossover scaling function beyond n = oo where the exponents <p/v and 7/V differ has yet to be 
investigated. 
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FIG. 1. Semi- logarithmic plot of the two-dimensional correlation length, showing the scaling as a function of temperature 
in the case of nesting. 

Fig. 2 
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FIG. 2. Neel Temperature as a function of t z = t\\ for U = 4t±, at half-filling. 
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Fig. 3 
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FIG. 3. as a function of ^ N |ln ^-| for U — 4t± at half-filling. The quantities T and t z = t\\ are in units of t± — t. 




FIG. 4. In-plane correlation length £± (lattice space is unity) for several values of out-of-plane hopping parameter at half 
filling for U = 4t±. T and t z = t\\ are in units of t±. 
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FIG. 5. Second in-plane derivative of the non-interacting susceptibility / (n, T) = d 2 xo (q,0) /d 2 q±\ci in two-dimensions as 
a function of temperature for various band-fillings n. 



Fig. 6 




FIG. 6. Plot of th e uni versa l crossover function from d = 2 to d = 3 for the staggered susceptibility in the n — > oo limit as 
defined by Eqs.(|D22|) and (|D20|). 
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